Optical Superradiance from Nuclear Spin Environment of Single Photon Emitters 
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We show that superradiant optical emission can be observed from the polarized nuclear spin 
ensemble surrounding a single photon emitter such as a single quantum dot (QD) or Nitrogen- 
Vacancy (NV) center. The superradiant light is emitted under optical pumping conditions and 
would be observable with realistic experimental parameters. 
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Superradiance (SR) is a cooperative radiation effect re- 
sulting from spontaneous build-up and reinforcement of 
correlations between initially independent dipoles. Its 
most prominent feature is an emission intensity burst 
in which the system radiates much faster than an oth- 
erwise identical system of independent emitters. This 
phenomenon is of fundamental importance in quantum 
optics and since its first prediction by Dicke in 1954 [l[ it 
has been studied extensively (for a review see 0)- The 
rich steady state properties of the associated dynamics 
can account for strong correlation effects including phase 
transitions and bistability 0, 0] ■ Yet in its original form 
optical SR is difficult to observe due to dephasing dipole- 
dipole Van der Waals interactions, which suppress a co- 
herence build-up in atomic ensembles. 

In this Letter we show that cooperative emission can 
occur from the ensemble of nuclear spins surrounding a 
quantum emitter such as a self-assembled QD or an NV 
center. The interaction of nuclear spin ensemble and op- 
tical field is mediated by the electron spin of the emitter. 
Due to the indirect character of the interaction, the de- 
phasing Van der Waals interactions vanish in this setting. 

We first explain the proposal using the example of an 
NV center in diamond. Then, we adapt the model to 
QDs, which promise strong effects due to the large num- 
ber of involved nuclei. Despite the inhomogeneity of the 
nuclear spin coupling and related dephasing processes, 
we predict a SR-like correlation build-up in the nuclear 
spin ensemble and a significant intensity burst of sev- 
eral orders of magnitude in the optical emission profile. 
Finally, we point out the possibility of observing phase 
transitions and bistability in the nuclear system. 

The superradiant effect is based on the collective hy- 
pcrfine (HF) interaction of the electronic spin of the de- 
fect (QD or NV) with N initially polarized proximal nu- 
clear spins. It is dominated by the isotropic contact term 
0, Q and reads in an external magnetic field (H — 1): 



H 



-(A+S- + A-S+) + gA z S z + u s S z 



(1) 



Here and A^ — J2i=i9i a i (/•* = +i ~ > z ) denote 
electron and collective nuclear spin operators, respec- 
tively. The coupling coefficients are normalized such that 



J2i 9i = 1 and individual nuclear spin operators are 
assumed to be spin-1/2 for simplicity; g gives the over- 
all HF coupling strength and wg denotes the electron 
Zeeman splitting. We neglect the typically very small 
nuclear Zeeman and nuclear dipole-dipole terms. 

Let us first consider NV centers, in which the effect 
can be studied in a clean and relatively small spin en- 
vironment. Due to their extraordinary quantum proper- 
ties, such as ultra-long decoherence times even at room 
temperature, NV centers have attracted wide interest Q 
resulting, e.g., in the demonstration of entanglement and 
quantum gates between the electron and proximal nu- 
clear spins ||. Both the NV center's electronic ground 
( 3 ^4) and optically excited states ( 3 E) are spin triplet 
(5=1) 7]. In the absence of a magnetic field, the ground 
state sublevels \ms = ±1) are split from \ms = 0). In 
the following, we assume that a B field is applied along 
the NV axis to bring \ms — 0) and \ms = 1) close to 
degeneracy |9J. In this case \ms = — 1) is off-resonance 
and can be disregarded. We focus on low-strain NV cen- 
ters with well-defined selection rules and assume that it 
is optically excited by selectively driving the weakly al- 
lowed transition from \ms = 1) to a state \E X ) in the 3 E 
manifold which decays primarily into |mg = 0) [Io| . see 
Fig. DJa) . The nuclear spin environment of the NV center 
consists of proximal 13 C (/ = 1/2) nuclei in the other- 
wise spinless 12 C matrix, which are HF-coupled to the 
electronic spin of the defect center. The interaction is 
dominated by the Fermi-contact term such that the cou- 
pling is isotropic (to first order) and described by Eq. ^ 
(ojs here contains both zero field splitting and Zeeman 
energy). Nevertheless in the simulations conducted below 
we included the small anisotropic dipole-dipole terms. 

We describe now the coupling of the nuclear spin to 
the optical field as depicted in Fig. [ljb) . It is best un- 
derstood as a two-step process: first, strongly driving a 
dipole- forbidden optical transition of the \ms = 1) spin 
state (the allowed transition is far off-resonant) pumps 
the electron into \ms = 0). Such Raman spin-flip tran- 
sitions have been demonstrated recently [10|. Since the 
short-lived excited state is populated negligibly through- 
out the process, we can eliminate it from the dynamics 
using standard techniques and obtain a master equation 
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for the electron spin decaying with effective rate T r 
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p = T r (S-pS + - -S + S-p - -pS + S-) - i[H, p], (2) 

each decay being accompanied by a Raman photon. Sec- 
ond, the return to state \ms = 1), necessary for the next 
emission, occurs through H via a HF mediated electron 
spin flip (and nuclear spin flop). Thus, each Raman pho- 
ton indicates a nuclear spin flop and the emission in- 
tensity I(t) is proportional to the change in nuclear po- 
larization. Starting from a fully polarized state, SR is 
due to the increase in the operative HF matrix element 
(A + A~). The scale of the coupling is set by A := gJ2i 9i- 
For a small relative coupling strength e = A/ (2 A), where 
A := |r. r /2 + iuis\, the electron is predominantly in its 
\ m s — 0) spin state and we can project Eq. ([2]) to the re- 
spective subspace. The reduced master equation for the 
nuclear density operator reads 



p=Cr(A P A + -^A + A p 



\pA + A-) 



(3) 



■id[A + A ,p] -igm s [A z ,p], 



where c r = g 2 /{2A) 2 Y r and a = g 2 /(2A) 2 uj s . 

As the electron is optically pumped into \ms = 0), the 
last term - representing the electron's Knight field - in 
Eq. ([3]) vanishes. Assuming resonance (wg = 0) the 
equation closely resembles the SR master equation which 
has been discussed extensively in the context of atomic 
physics [3] and thus SR effects might be expected. How- 
ever, there is a crucial difference: the inhomogeneous 
nature (gi ^ const) of the operators A^: they do not 
preserve the collective spin, affecting the relative phase 
between nuclei. This could prevent the phased emission 
necessary for SR [1, HH, [Hj] . However, as we shall see, SR 
is still clearly present in realistic inhomogeneous systems. 
We take the ratio of the maximum intensity to the initial 
intensity (the maximum for independent spins) I C oop/^ind 
as our figure of merit in the following: if this relative in- 
tensity peak height is > 1 it indicates cooperative effects. 

To see that this effect can be observed at NV centers, 
we simulate Eq. ([2]) numerically jl3j . The number N of 
effectively coupled nuclei can range from a few to a few 
hundred, since the concentration of 13 C can be widely 
tuned [la]. The HF constants gi between the defect and 
the nearest ~ 40 nuclei were derived in Q in an ab-initio 
calculation. Nuclei outside this shell (~ 7 A) have a cou- 
pling strength ggi weaker than 2ix ■ 0.5MHz and are not 
considered here. The excited state lifetime of the NV 
center has been measured as r » 13ns 1(3, 17 1. Thus, we 
adopt an effective rate Y r = 2-k ■ 10MHz for the decay 
from \ms = 1) to \ms = 0) enabled by driving the Ra- 
man transition. The intensity enhancements predicted 
by exact simulations for small, randomly chosen and ini- 
tially polarized spin environments are shown in Fig.JIJc). 
In samples of higher 13 C concentration N can be larger 
and stronger effects are expected. 
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FIG. 1. (a) Simplified level scheme of NV center with relevant 
A system (cf. text and [10(); (b) Sketch of relevant processes: 
electronic ground states are coupled by optical pumping and 
HF flipflops; the states are labeled by the ^-components of the 
electron and nuclear spin, (c) 7 C oo P //ind for randomly chosen 
nuclear environments of an NV center The first nuclear shell 
is taken to be spinless, as due to their very strong coupling 
they would evolve largely independent from the ensemble. 



One characteristic feature of SR is the linear TV- 
dependence of the associated effects (already visible in 
Fig- Hi- Since the number of nuclei to which the electron 
couples is much larger in a QD than in a NV center, QDs 
are particularly attractive candidates for the investiga- 
tion of SR. In the following we study the dynamics of the 
QD system in different regimes and we show that strong 
signatures of SR can be expected in realistic settings. 

Let us consider a self-assembled QD in which a single 
conduction band electron is coupled by isotropic Fermi 
contact interaction to a large number of nuclear spins. 
Optical pumping of the electron is realized by a Raman 
process, driving a forbidden transition to a trion state 
[18j |. and including the HF coupling we again obtain dy- 
namics as sketched in Fig.QJb). For the optical pum pin g 
rate values T r = 2tt ■ (0.1 — l)GHz are applicable [l9ll20|. 
A comparison with the HF coupling constants reported 
for different materials [l4[ shows that for InGaAs and 
CdSe QDs at resonance Eq. ([3]) is not valid since the 
relative coupling strength e > 1. We therefore consider 
the dynamics of the system under conditions of a finite 
electron inversion [using Eq. ((5)]. In this regime, the 
electron can be seen as a driven and damped two-level 
system: the nuclei 'pump' excitations into the electron, 
which are damped by the Raman-mediated decay; coop- 
erative behavior manifests in enhanced HF interaction. 
This enhancement directly translates into increased elec- 
tron inversion (S + S~) to which the emitted photon rate 
is proportional and thus SR from a single QD can be 
expected. Let us rephrase this, since SR from a single 
emitter is somewhat counter-intuitive. Of course, on an 
optical time scale, anti-bunched single photons will be 
emitted at a rate below the optical decay rate. It is, in 
fact, typically much slower since the emitter is pumped 
into the optically inactive state \ms = 0). SR on time 
scales ~ 1 /r r consists thus of lifting this "spin-blockade" 
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by HF coupling which becomes increasingly more efficient 
as nuclear cooperative effects kick in. As in the homo- 
geneous case [2J this enhancement is associated with the 
transition through nuclear Dicke states \J, m) , \m\ <C J. 
Although J is not preserved by inhomogeneous we 
can use the Dicke states to illustrate the dynamics. For 
instance, due to the large homogeneous component in 
A~ , its matrix elements show a strong increase oc J for 
states | J, m ) , |m| <C J. 

For large relative coupling strengths e>l the electron 
saturates and superradiant emission is capped by the de- 
cay rate T r /2, prohibiting the observation of an intensity 
burst. In order to avoid this bottleneck regime we choose 
a detuning uj s = A/2 such that < e = A/y/T 2 . + A 2 < 
1. In this parameter range the early stage of the evolution 
- in which the correlation build-up necessary for SR takes 
place Q - is well described by Eq. ([5} . The nuclear phas- 
ing is counteracted by the dephasing (inhomogeneous) 
part of the Knight term (oc g^/Var(<?i)/2 (HI), which can 
cause transitions J — > J — 1. However, the system evolves 
in a many-body protected manifold (MPM) 122]: The 
term ~ [A + A~,p] energetically separates different total 
nuclear spin- J manifolds. A rough estimate of the ra- 
tio between detuning and dephasing shows a dependence 
oc e 2 , with proportionality factor > 1 (diverging in the 
homogeneous limit). Thus for values e « 1 the correla- 
tion build-up should be largely MPM-protected. We now 
confirm these considerations and show by numerical sim- 
ulation of Eq. ([2]) that a SR peaking of several orders of 
magnitude can be observed in the Raman radiation from 
an optically pumped QD, cf. Fig. [2j An exact numeri- 
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FIG. 2. Relative intensity under dynamical Overhauser field 
compensation: N = 21 2 , lj s = A/2 and e = 0.3(1), 0.7(11), 
0.99(111). (IV) shows the ideal Dicke SR profile [U as a refer- 
ence. Inset: comparison of exact evolution and factorization 
for TV = 9 inhomogeneously coupled spins (left peak, scaled 
by factor 50) and N = 21 2 homogeneous spins (e = 0.7). Fully 
independent emitters lead to an exponential curve slowly de- 
caying from 1 to zero and are therefore not depicted. 



cal simulation of the dynamics is not feasible due to the 
large number of coupled nuclei and since the dynamics for 
inhomogeneous coupling cannot be restricted to a low- 
dimensional subspace. To obtain I(t) oc 



we therefore use an approximative scheme. By Eq. ([2]) 
these expectation values are related to fourth-order cor- 
relation terms involving both the electron and nuclear 
spins. We use a factorization assumption to reduce the 
higher-order expressions in terms of the covariance ma- 
trix 7^ = ^cr+CT^). Following [23[ we apply the bosonic 
Wick's theorem, incorporating the fermionic character of 
same-site nuclear spin operators ([<7T , = 1) an d rc_ 

place, e.g., (a+a^-) ( 7 + - h)(a+S~) 



e.g., {a i ap ) -> (j^ - b ) - 74 x _ , 

But the electron spin plays a special role and factoriz- 
ing it completely leads to poor results. Therefore we 
also solve Eq. ([5]) for the main higher-order term in- 
volving the electron, the "mediated covariance matrix" 
7,^ = (cr^S z o-j). All other higher-order expectation val- 
ues therein are factorized under consideration of special 
symmetries for operators acting on the same site. 

In the regimes accessible to an exact treatment, i.e., 
the homogeneous case and for few inhomogeneously cou- 
pled particles, the factorization results agree well with 
the exact evolution (see inset in Fig. [2]). This shows that 
it quantitatively captures the effect of nuclear spin co- 
herences while allowing a numerical treatment of hun- 
dreds of spins. Finally, in addition to the constant de- 
tuning us = A/2 for the displayed simulations we com- 
pensated the Overhauser field dynamically 24[ . Further- 
more, we assume a Gaussian spatial electron wave func- 
tion. The results obtained with these methods are dis- 
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FIG. 3. /coop/^ind for different values of e - the Overhauser 
field is dynamically compensated and ujs = A/2 in all cases 
- compared to the ideal Dicke case. Dashed line corresponds 
to a partially polarized dot (cf. text). 



played in Figs. [2] and [3] For e w 1 the strong MPM 
protection suppresses dephasing, leading to pronounced 
SR signatures: A strong intensity burst, whose relative 
height scales oc N (for large N). For large e w 1 the 
relative height is reduced only by half compared to the 
ideal Dicke case. For smaller e it decreases further due 
to increased dephasing. For e < 0.3, where MPM protec- 
tion is weak and the decay process is significantly slowed 
down (c r oc e 2 ), even the linear scaling is lost. 

From Fig. [3] one extrapolates that for a fully polarized 
initial state a huge intensity overhead of several orders of 



4 



magnitude (~ 10 3 -10 4 ) is predicted. If the initial state 
is not fully polarized, SR effects are reduced. However, 
even when, e.g,. starting from a mixture of symmetric 
Dicke states | J, J) with polarization P = 60% [3,[2^] our 
simulations predict a strong intensity peak and (for TV 3> 
1) a linear TV-dependence: I COO p/hnd ~ 0.037V (e = 0.99), 
i.e., only a factor 4 weaker than for full polarization. 

Note that for the sake of simplicity we consider 1=1/2 
nuclei in our simulations. In terms of particle numbers 
TV this is a pessimistic assumption as typical QD host 
materials carry a higher spin. We can incorporate this 
effect by treating higher spins as 21 homogeneously cou- 
pled spins 1 = 1/2 thus increasing the effective particle 
number by the factor 21. Most QDs consist of a few dif- 
ferent species of nuclei with strongly varying magnetic 
moments, increasing the inhomogeneity of the system. 
However, in the worst case the different species evolve 
independently diminishing the effect by a small factor 
corresponding to the number of species. In our simula- 
tions the effect was shown to be much smaller. 

We have neglected the dipolar and quadrupolar inter- 
action among the nuclear spins. The former is always 
negligible on the time scale considered here [5|. The lat- 
ter is absent for nuclear spin 1 = 1/2 (NV centers, CdSe 
QDs) or strain-free QD s [19j ]. In strongly strained QDs 
it can be important [261 ] . and a term ^(Z?*) 2 must be 
added to Eq. ([1]), where z% is defined as the main axis of 
the local electric-field-gradient tensor. 

Having seen that SR can be observed in experimentally 
accessible nuclear spin ensembles, let us briefly explore 
two further applications of this setting: nuclear spin po- 
larization and phase transitions. We first note that the 
master equation Eq. ^ describes optical pumping of the 
nuclear spins. Its steady states are the eigenstates of A z 
which lie in the kernel of A~ , so-called dark states, and 
include the fully polarized state. Hence the setting de- 
scribed by Eqs. (|2|3p can be used to polarize the nuclei 
[13], i.e., to prepare an initial state required for SR. 

Finally, nuclear spin systems may be used to study fur- 
ther cooperative effects such as phase transitions. It is 
known [3j that in the thermodynamic limit an optically 
driven atomic system with collective decay - as described 
by Eq. (J3j> for homogeneous operators and ms = = 
- can undergo a second-order non-equilibrium phase 
transition in the steady state. In our setting an effec- 
tive driving can be established by a DC magnetic field 
B x perpendicular to the polarization direction. A semi- 
classical treatment of the equations of motion deduced 
from Eq. (|2|) predicts a similar phase transition in the 
combined system of electron and nuclear spins in certain 
regimes. Preliminary simulations confirm the validity of 
the semiclassical results and also indicate the appearance 
of related phenomena like bistability and hysteresis which 
have recently been observed in polarization experiments, 
e.g. [m . A detailed analysis of these topics and an an- 
alytical description of the SR dynamics presented here 



will be subject of a forthcoming publication [29(. 

In conclusion, we have shown that the nuclear spin en- 
vironment of individual QDs and NV centers shows su- 
perradiant optical emission under suitable optical pump- 
ing conditions. While in NV centers a collective intensity 
enhancement of up to 100% is predicted, the much larger 
nuclear spin ensembles in QDs could lead to relative peak 
heights of several orders of magnitude. This would be 
clear evidence of coherent HF dynamics of nuclear spin 
ensembles in QDs. The rich physics of SR systems, in- 
cluding bistability and phase transitions could thus be 
studied in a long-lived mesoscopic solid-state system. 
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